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INTRODUCTION 


Previous  communications  from  this  laboratory  have  reported  results 
of  calculations  on  the  adiabatic  flame  temperatures  and  resulting  pres¬ 
sures  from  the  constant-volume  combustions  of  organic  fuels  in  air^  and 
of  magnesium  in  air.^  The  present  study  represents  a  continuation  of 
the  foregoing,  considering  the  effects  to  be  observed  when  magnesium 
and  fuel  together  undergo  adiabatic  constant-volume  combustion  in  air. 

As  was  the  case  previously,  the  calculations  have  been  carried  out 
on  a  desk-top  computer  (WANG  720C) ,  Magnesium  and  fuel  are  considered 
to  be  uniformly  distributed  in  the  air  volume  at  25 °C,  one  atmosphere. 

The  air  is  assumed  to  be  78  mole  percent  N2,  21  mole  percent  O2,  and 
One  mole  percent  Ar. 

The  fuels  considered  were  mannitol  hexanitrate  (nltromannite) ,  PETN, 
HMX,  Composition  B-3,  TNT,  ethylene  oxide,  and  THDC  (exo-tetrahydrodi- 
cyclopentadiene) .  Formulas  and  a  few  significant  properties  of  the  fuels 
are  given  in  Table  1. 

In  all  the  following,  the  symbol  C  will  represent  the  charge,  or 
the  mass  of  fuel  per  cubic  meter  of  air,  and  M  the  mass  of  magnesium  per 
cubic  meter  of  air.  Charge-to-mass  (C/M)  ratios  of  10:1,  3:1,  1:1,  1:3, 
and  1:10  were  examined,  with  total  concentration,  C  +  M,  ranging  from 
0.02  kg/m^  to  10  kg/m^. 


^  Naval  Weapons  Center.  Peak  Overpressures  for  Internal  Blastyhj 
G.  F.  Kinney,  R.  G.  S.  Sewell,  and  K.  I.  Graham.  China  Lake,  Calif.', 
NWC.  (NWC  TP  6089,  publication  UNCLASSIFIED.)  (In  press.) 

^  Naval  Weapons  Center.  Reaoti-ve  Metals  in  Internal  Explosions: 
The  Combustion  of  Magnesium  in  Air^  by  R.  A.  Reinhardt.  China  Lake,. 
Calif.,  NWC,  February  1978.  (NWC  TM  3429,  publication  UNCLASSIFIED.) 
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TABLE  1.  Properties  of  the  Fuels 


Designa¬ 

tion 

Name  and  Chemical  Formula 

^^f ,298 
(kJ/mole) 

Relative  Oxygen 

Content^ 

MHN 

Mannitol  hexanitrate  (nitromannite) 

-608.96 

1.80 

PETN 

Pentaerythritol  tetranitrate . 

SV4"l2 

-489.8 

1.33 

HMX 

sym-Cyclotetramethylenetetra- 
nitramine ,  C^HgNgOg 

+77.4 

1 

1.00 

Comp  B— 3 

60%  RDX^,  40%  TNT 

r  M  w  n  ^  ! 

'"2*05  2.51^2.15  2.67  I 

-0.82'^ 

0.81 

TNT 

2,4,6-Trinitrotoluene 

SS^3°6 

-57.25 

0.63 

Ethylene 

Oxide 

Epoxyethane 

C^H^O 

-47.29 

0.25 

THDC 

exo-Tetrahydrodicyclopentadiene 

^10\6 

-104.9 

0.00 

^  Number  of  oxygen  atoms/ (number  of  carbon  plus  half  the  number 
of  hydrogen  atoms). 

sym~Trlmethylenetrinitramine  (cyclonite) . 

^  For  a  nominal  formula  mass  of  100  g/mole. 
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BASIS  OF  CALCULATIONS 


From  the  First  Law  of  Thermodynamics  it  follows  that  for  an  adia¬ 
batic,  constant-volume  non-flow  process  there  must  be  no  net  change  in 
U,  the  internal  energy;  or  the  internal  energy  of  all  the  products  at 
the  adiabatic  flame  temperature,  T,  must  equal  that  of  the  reactants  at 
entering  conditions.  For  each  component  in  the  mixture  of  products, 
using  the  ideal-gas  assumption  throughout,  we  may  write: 

T 

Ut  “  1^298  ^ 

298 

where  Cy  is  the  heat  capacity  at  constant  volume.  For  U298  the  custom¬ 
ary  thermochemical  convention  may  be  chosen  of  using  the  internal  energy 
of  formation  of  298  K.  The  latter  is  defined  as  the  internal  energy 
change  for  the  synthesis  of  the  substance  considered  from  the  elements 
in  their  stable  states  at  the  given  temperature  (298  K) .  Of  the  re¬ 
actant  materials  all  but  the  organic  fuel  are  elements  in  the  stable 
state  at  298  K,  and  thus  the  initial  Internal  energy  is  just  equal  to 
the  Internal  energy  of  formation  of  the  fuel,  AUf,  at  298.  Hence,  the 
adiabatic-isochoric  condition  may  be  written: 

E  Ui,T  = 


the  index,  i,  applying  to  all  the  product  species. 

A  temperature  must  be  found  for  which  this  criterion  is  met  and 
which  may  then  be  identified  as  the  adiabatic  flame  temperature.  The 
pressure  in  the  system  is  then  calculated  from  the  ideal  gas  law: 


P  =  RT/V  ^n^ 

in  which  n^  is  the  mole  number  for  the  ith  product  component ,  consider¬ 
ing  only  the  gaseous  products,  V  the  total  volume  (one  cubic  meter)  and 
R  the  universal  gas  constant .  The  overpressure  is  found  by  subtracting 
one  bar. 

As  before,^  thermochemical  data  are  represented  by  a  five-parameter 
expression: 

Ut  =  A  +  BT  +  CT^  +  DT^  +  E  log  T 
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Table  2  gives  the  internal  energy  parameters  for  the  products  considered- 
data  from  the  JANAF  Tables 3  were  used  as  the  primary  source. 

In  Appendix  C  is  described  the  condensed  storage  of  the  parameters 
pressing  equilibrium  constant  and  internal  energy  as  functions  of  tern- 

psrs  tTpUTB « 

equilibrium  considerations 

The  following  species  are  considered  present  in  the  product  mixture: 
solid  or  liquid  MgO  (depending  on  whether  the  temperature  is  above  or 
below  the  melting  point,  3050  K) ,  perhaps  solid  carbon,  and  the  thirteen 
gases  CO,  CO2,  H,  OH,  H2O,  t(gO,  0,  NO,  Ar,  H2,  Mg,  N2,  and  O2.  Speaking 
generally,  in  order  to  solve  for  the  fifteen  mole  numbers  of  these 
species,  there  must  be  fifteen  independent  relations.  Of  these,  6  are 
the  atom  balance  conditions  (see  Appendix  A)  for  the  6  elements  repre¬ 
sented,  and,  thus,  9  equilibrium  relations  are  required.  These  may  be 
considered  to  be  the  vapor  pressures  of  MgO  and  the  equilibrium  con- 

Sefnara  “2*  H,  OH,  H2O,  MgO(g),  0.  and  NO.  The 

three-parameter  equation  used  by  Kinney,  et  alj^log  K  =  A  +  B/(C  +  T) 

based  on  data  from  the  JANAF  Thermochemical  Tables 3  has  been  used.  In 

values  of  the  parameters  A,  B,  and  C  which  were 

waf  a  t  f  .The  methods  for  determining  these  parameters 

was  described  previously,^ 

f  a  computation,  at  an  assumed  temperature,  T,  is  greatly  simpli¬ 
fied  by  making  use  of  the  chemical  behavior  of  the  system.  The  follow¬ 
ing  outlines  the  method  used.  More  detail  is  given  in  Appendix  B. 

1  the  vapor  pressure  of  the 

solid  (if  below  3050  K,  the  melting  point)  or  liquid  at  the  temperature 

h-  conclusion  of  computation  showed  no  case  in 

which  all  the  MgO  present  was  vaporized.  In  the  most  general  case, 
there  are  6  gaseous  species  to  consider  first  in  the  computation  (the 
or  species _  of  Appendix  B) :  H2,  H2O,  CO,  CO2,  O2,  Mg.  From  an 
nitial  approximation  of  one  of  them,  the  remaining  5  may  be  found  from 
M  appropriate  combination  of  atom  balance  and  equilibrium  relations. 

The  resulting  5  simultaneous  nonlinear  equations  lead  to  a  cubic  equa¬ 
tion  in  one  unknown  which  may  be  solved  iteratively.  The  results  of 

computation  of  the  "high-temperature  spe¬ 
cies  ,  M,  U,  OH,  and  NO.  An  adjusted  atom-balance  conditions  leads  to 
the  next  iteration. 


3  National  Bureau  of  Standards.  JANAF  Thermo ahemioal  Tables  by 
Ed.^'june^Wyi*^  Directors.  Washington,  D.  C."  2nd 
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TABLE  2.  Internal  Energies  of  Products  of  Combustion 
of  Magnesium  plus  Organic  Fuel  in  Air 

Expressed  as  a  function  of  x,  the  absolute  temperature  in  kilo- 
kelvins  (x  =  T/1000) ;  values  given  in  joules/mole: 

2  3 

=  A  +  Bx  +  Cx  +  Dx  +  E  In  X 


Substance 

A 

B 

C 

D 

E 

CO 

-125179 

28870 

433 

-37.6 

-4560 

CO2 

-421544 

56400 

251 

13.4 

-10000 

H 

213050 

12470 

0 

0 

0 

OH 

33024 

18623 

2533 

-192.8 

-450 

H2O 

-254588 

29480 

5084 

-403.2 

-4930 

MgO(s)® 

-621030 

53360 

474 

117.8 

-3380 

HgOCl)"^ 

-562034 

60670 

0 

1 

0 

0 

MgO(g)  1 

-7635 

28690 

299 

-3.5 

0 

NO 

75341 

31460 

-118  ! 

9.6  1 

-5560 

0 

245728 

11700 

Ill 

8.8 

0 

Ar 

-3718 

12470 

0 

0 

0 

«2 

-3789 

15900 

2933 

-184.2 

870 

Mg(g) 

136547 

15270 

-2339 

326.9 

7980 

^2 

-16128 

31900 

-199 

15.9 

-7220 

O2 

-8976 

24270 

1671 

-88.3 

-770 

C(s) 

-17007 

29350 

-524 

52.9 

-6800 

^  Values  for  solid  and  liquid  MgO  have  been  slightly 
corrected  from  those  in  Reference  2. 
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TABLE  3.  Equilibrium  Constants  Involved  in 
Magnesium-Fuel-Air  Interactions 

Logarithms  (base  10)  of  the  equilibrium  constant  of  the  indicated 
reactions  expressed  as  a  function  of  absolute  temperature 

logj^Q  K  =  A  +  B/{C+T) 


Constant 

Reaction 

A 

B 

r - — 

c 

^CO 

C(s)  +  hO^  =  CO 

4.069 

7440 

184 

^C02 

C(s)  +  0^  =  CO2 

-0.099 

21100 

18 

Kr 

J5H2  =  H 

3.242 

-12280 

35 

%H2=  =  OH 

0.662 

-1720 

-122 

^2° 

h  +  *Stl2  -  HjO 

-3.100 

13420 

20 

OQ 

O 

Mg(g)  +  Ho^  =  MgO(g) 

-1.772 

-181 

^0 

JsN2  +  ho^  =  NO 

0.645 

im 

-15 

^0 

%02  =  0 

3.545 

-13550 

16 

MgO(s)  =  MgO(g) 

8.239 

-27570 

-104 

MgO(l)  =  MgO(g) 

6.65 

-22500 

-154 
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Several  special  cases  should  be  indicated: 

Case  I.  Oxygen  is  in  excess  over  both  carbon  (to  CO)  and  magnesium. 
Case  II.  Insufficient  oxygen  for  both  carbon  and  magnesium;  no  solid 
carbon  present . 

Case  III.  Solid  carbon  is  present. 

Case  IV.  At  the  extreme  of  high  magnesium  concentration,  crystalline 
magnesium  nitride  will  form,  sometimes  accompanied  by  condensed  ele¬ 
mental  magnesium  metal  (Table  4).  Because  of  the  memory  limitations^ 
of  the  WANG  calculator,  it  was  not  practicable  to  include  this  case  in 
the  computations.  These  points  have  thus  been  omitted  (and  are  so  in¬ 
dicated  in  Table  5) . 


RESULTS 

In  Table  5  are  given  the  ][;esults  of  the  computations.  The  table 
is  in  five  sections,  corresponding  to  the  five  values  of  the  charge-to- 
mass  ratio  (C/M)  chosen.  For  each  point  the  upper  figure  is  the  over¬ 
pressure  (P  -  1,  in  bars)  and  the  lower,  in  parentheses,  is  the  Kelvin 
temperature.  Moving  down  any  column  of  this  table  corresponds  to  in¬ 
creasing  the  total  concentration,  C  +  M,  for  a  fuel-magnesium  mixture 
of  fixed  C/M.  Typically,  the  temperature  reaches  a  maximum  about  where 
the  system  is  just  balanced  in  oxygen  for  the  formation  of  CO  +  H2O. 

For  example,  for  TNT  at  C/M  =  1.0,  the  temperature  maximum  occurs  at 
C  =  0.5  kg/m^,  M  =  0.5  kg/m^.  Beyond  here,  it  is  to  be  noted  that  the 
pressure  continues  to  rise  even  as  the  temperature  falls  as  a  result  of 
increasing  quantity  of  gaseous  products.  At  the  last  point  (C  and  M^ 
each  5  kg/m^)  the  temperature  is  a  bit  higher  than  the  preceding;  this 
is  always  observed  when  the  smoke  point  is  reached  (when  solid  carbon 
forms).  The  carbon  is  formed  chiefly  from  the  reaction  C0(g)  =  C(s)  + 
MgO(s),  a  highly  exothermic  reaction. 

The  points  indicated  by  note  ^  in  Table  5  are  those  for  which 
Mg3N2  is  produced.  As  pointed  out  earlier,  it  was  not  practicable  to 
program  for  this  case;  it  can  be  said,  nonetheless,  that  the  pressures 
will  be  very  much  smaller  than  otherwise  because  of  the  consumption  of 
four  moles  of  gas  per  mole  of  Mg3N2  formed. 

Table  6  gives  product  yields  (as  mole-fraction  in  the  gas)  for  this 
series  just  referred  to:  TNT  at  C/M  =  1.0.  The  extreme  complexity  found 
at  the  middle  concentrations— say ,  from  C  +  M=  0.20  to  1.0  kg/m  — 
should  be  noted:  each  of  the  thirteen  gaseous  species  is  present  in 
appreciable  quantity.  Such  data,  though  not  included  in  this  report, 
are  available  for  all  the  other  points  computed. 

Figures  1  and  2  show  the  pronounced  effect  of  total  concentration 
(C  +  M)  on  overpressure;  the  cases  chosen  for  Figure  1  are  for  HMX  at 
three  charge-to-mass  ratios.  About  the  same  appearance  is  to  be  ob¬ 
served  with  any  of  the  other  explosive  fuels.  The  hydrocarbon,  THDC, 
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TABLE  4.  Logarithm  (base  10)  of  the  Dissociation 
Constant  for  Magnesium  Nitride 


K 

P 


P 

Mg  N 


2 


T,  K 

^°®10 

900 

-19.56 

1000 

-16.21 

1100 

-11.10 

1200 

-9.16 

1300 

-7.52 

1400 

-5.89 

1500 

-3.71 

1600 

-1.82 

1700 

-0.15 

1800 

1.33 

1900 

2.65 

2000 

3.84 

2100 

4.90 

2200 

5.87 

2300 

6.75 
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TABLE  5.  Overpressure  (bars)  and  Temperature  (in  parentheses) 
for  the  Adiabatic,  Constant -Volume  Combustion  of  Fuels 
Plus  Magnesium  in  Air. 


C/M 

Fuel 

r 

Mg 

NM 

PETN 

HMX 

B-3 

TNT 

^2\° 

THDC 

l:l 

.01 

.01 

1.16 

(638) 

1.22 

(659) 

1.26 

(671) 

1.34 

(695) 

1.46 

(732) 

1.92 

(866) 

2.39 

(1006) 

.02 

.02 

2.19 

(945) 

2.29 

(977) 

2.39 

(998) 

2.51 

(1040) 

2.71 

(1103) 

3.49 

(1329) 

4.28 

(1566) 

.05 

.05 

4.7 

(1670) 

4.9 

(1732) 

5.0 

(1774) 

5.3 

(1856) 

5.7 

(1980) 

7.1 

(2369) 

8.3 

(2700) 

.1 

.1 

7.7 

(2526) 

8.0 

(2601) 

8.3 

(2657) 

8.5 

(2752) 

8.9 

(2869) 

10.1 

(3050) 

10.7 

(3034) 

.2 

.2 

10.4 

(3107) 

10.7 

(3147) 

11.0 

(3193) 

11.2 

(3234) 

11.5 

(3276) 

12.9 

(3149) 

10.3 

(2183) 

.5 

.5 

16.2 

(3574) 

16.8 

(3559) 

17.7 

(3550) 

17.8 

(3503) 

17.9 

(3371) 

16.5 

(2297) 

13.7 

(2201) 

1 

1 

25.4 

(3766) 

26.5 

(3699) 

28.2 

(3618) 

28.0 

(3457) 

27.0 

(3050) 

25.1 

(2363) 

18.2 

(2010) 

2 

2 

44.0 

(3922) 

46.1 

(3802) 

49.1 

(3636) 

47.2 

(3284) 

39.2 

(2501) 

42.4 

(2441) 

a 

5 

5 

101.1 

(4138) 

105.9 

(3948) 

112.1 

(3652) 

103.2 

(3050) 

87.7 

(2655) 

95.3 

(2566) 

a 

1:3 

.005 

.015 

1.47 

(735) 

1.51 

(750) 

1.53 

(756) 

1.56 

1  (767) 

1.61 

(785) 

1.85 

(852) 

2.09 

(923) 

.01 

.03 

2.72 

(1118) 

2.77 

(1133) 

2.81 

(1143) 

j  2.87 

1  (1163) 

2.97 

(1194) 

3.35 

(1306) 

3.75 

(1428) 

.02 

.06 

4.8 

(1742) 

4.8 

(1767) 

4.9 

(1784) 

j  5.0 
(1818) 

5.2 

(1870) 

1 

1 

J - 

5.8 

(2046) 

- 

6.4 

(2236) 
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TABLE  5.  (Cont'd) 


C/M 

Fuel 

Mg 

1 - 

NM 

PETN 

HMX 

B-3 

TNT 

C2H4O 

THDC 

1:3 

.05 

.15 

8.6 

(2949) 

8.7 

(2967) 

8.8 

(2988) 

8.9 

(3016) 

9.0 

(3050) 

9.2 

(3050) 

9.7 

(3128) 

.1 

.3 

10.6 

(3378) 

10.7 

(3382) 

10.9 

(3394) 

11.0 

(3402) 

11.1 

(3411) 

12.0 

(3364) 

12.0 

(3109) 

.2 

.6 

.  14.3 
(3613) 

14.5 

(3588) 

14.9 

(3568) 

15.0 

(3538) 

15.1 

(3482) 

15.8 

(3132) 

12.7 

(2260) 

.5 

1.5 

25.2 

(3648) 

25.5 

(3557) 

26.0 

(3440) 

25.4 

(3261) 

25.3 

(3050) 

22.5 

(2363) 

19.0 

(2194) 

1 

3 

43.0 

(3588) 

42.8 

(3398) 

41.8 

(3087) 

41.3 

(2885) 

35.4 

(2461) 

36.8 

(2408) 

a 

2 

6 

77.7 

(3485) 

74.4 

(3138) 

76.2 

(2931) 

66.7 

(2551) 

63.2 

(2543) 

65.3 

(2439) 

a 

1:10 

.002 

.02 

' 

1.79 

(848) 

1.83 

1  (851) 

1.84 

(854)  i 

1 

1.85 

(858) 

1.87  : 

(865) 

1.96 

(891) 

2.06 

(919) 

.004 

.04 

3.29 

(1301) 

3.31 

(1306) 

3.32 

(1310) 

3.35 

(1318) 

3.38 

(1330) 

3.50 

(1374) 

3.69 

(1422) 

.01 

.1 

6.6 

(2352) 

6.6  ■ 
(2361) 

6.6 

(2369) 

6.7 

(2382) 

6.7 

(2405) 

7.0 

(2473) 

7.2 

(2552) 

.02 

.2 

8.6 

(3075) 

8.5 

(3077) 

8.7 

(3087) 

8.7 

(3097) 

8.8 

(3112) 

9.0 

(3144) 

9.2 

(3184) 

.04 

.4 

11.0 

(3514) 

11.0 

(3510) 

11.6 

(3521) 

11.1 

(3510) 

11.2 

(3509) 

11.5 

(3479) 

11.8 

(3424) 

.1 

1 

1 

1 

17.0 

(3543) 

17.0 

(3515) 

17.2 

(3486) 

17.1 

(3454) 

17.1 

(3404) 

17.2 

(3183) 

16.5 

(2740) 
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TABLE  5.  (Cont'd) 


C/M 

Fuel 

~T 

Mg 

NM 

PETN 

B 

B-3 

TNT 

"2^4° 

THDC 

.2 

2 

( 

25.3 

:3201) 

24.9 

(3100) 

25.2 

(3050) 

25.7 

(3050) 

25.1 

(2906) 

21.8 

(2335) 

20.2 

(2244) 

.4 

4 

1 

39.9 

(2729) 

37.4 

(2500) 

35,.  4 
(2336) 

34.8 

(2327) 

— 

34.1 

(2318) 

33,3 

(2179) 

a 

3:1 

,015 

.005 

0.79 

(530) 

0.90 

(561) 

0.98 

(580) 

1.10 

(620) 

1.28 

(675) 

1.98 

(880) 

2.68 

(1087) 

.03 

.01 

1.60 

(761) 

1.78 

(813) 

1.87 

(846) 

2.11 

(912) 

2,44 

(1010) 

3.63 

(1351) 

4,79 

(1699) 

.06 

.02 

2.96 

(1139) 

3.26 

(1225) 

3.47 
j  (1282) 

3.84 

(1394) 

4.38 

(1560) 

6.3 

(2100) 

8.0 

(2568) 

1  .15 

.05 

6.2 

(1978) 

6.8 

(2128) 

7.2 

(2230) 

7.8 

(2419) 

8.6 

(2656) 

10.7 

(2977) 

9.9 

(2527) 

.3 

.1 

10.1 

(2766) 

10.7 

(2885) 

11.3 

(2973) 

11.5 

(3050) 

11.8 

(3113) 

13.5 

(2897) 

9.9 

(2055) 

.6 

.2 

14,3 

(3170) 

15.2 

(3233) 

16.1 

(3302) 

16.4 

(3324) 

16,6 

(3266) 

15.6 

(2217) 

12.0 

(2033) 

1.5 

,5 

26,1 

(3547) 

27,..8 

(3553) 

30.2 

(3562) 

30.2 

(3451) 

28.8 

(3050) 

27.3 

(2299) 

17.5 

(1803) 

3 

1 

45.4 

(3762) 

48.8 

(3733) 

53.5 

(3690) 

52.6 

(3468) 

42.3 

(2481) 

47.0 

(2376) 

a 

6 

2 

84.6 

(3962) 

91.2 

(3904) 

100.7 

(3817) 

97.4 

(3474) 

77.9 

(2597) 

87.0 

(2466) 

a 
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TABLE  5.  (Cont'd) 


C/M 

1 

Fuel 

Mg 

1 - 

NM 

PETN 

■ 

HMX 

B-3 

TNT 

S»4° 

THDC 

10:1 

.02 

.002 

0.61 

(473) 

0.76 

(516) 

0.85 

(544) 

1.03 

(598) 

1.29 

(676) 

1.72 

(941) 

3.11 

(1208) 

.04 

.004 

1.30 

(667) 

1.55 

(739) 

1.71 

(784) 

2.02 

(874) 

2.45 

(1006) 

4.02 

(1451) 

5.52 

(1893) 

.10 

.01 

3.01 

(1114) 

3.51 

(1253) 

3.87 

(1344) 

4.46 

(0522) 

5.3 

(1782) 

8.2 

(2514) 

9.2 

(2713) 

.2 

.02 

5.4 

(1667) 

6.2 

(1881) 

6.9 

(2024) 

7.8 

(2228) 

8.9 

(2611) 

11.1 

(2873) 

8.2 

(1894) 

.4 

.04 

9.3 

i  (2379) 

10.4 

(2615) 

11.3 

(2763) 

12.0 

(2939) 

12.6 

(3050) 

13.0 

(2470) 

7.8 

(1507) 

1 

.1 

17.3 

(3050) 

18.5 

(3135) 

20.1 

(3234) 

20.5 

(3250) 

20.0 

(3067) 

17.8 

(2002) 

9.4 

(1226) 

2 

.2 

28.3 

(3291) 

30.7 

,(3363) 

33.8 

(3428) 

33.7 

(3331) 

31.2 

(2917) 

27.6 

(1945) 

12.3 

(1116)b 

4 

.4 

50.5 

(3494) 

55.6 

(3538) 

61.2 

(3571) 

59.8 

(3374) 

47.3 

(2434) 

47.2 

(1906) 

c 

10 

1 

118.3 

(3716) 

129.7 

(3738) 

144.6 

(3742) 

138.8 

(3440) 

109.0 

(2578) 

106.4 

(1884) 

a 

MgsN^  is  found  in  these  cases , 

1.8  kg/m  of  fuel,  0.18  k/m^  of  air. 
Failed  to  converge . 
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Mole  percent 


TABLE  6.  Product  Yield  for  TNT-Mg ,  C/M  -  1.0 
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FIGURE  2.  THDC:  Overpressure  vs  Total  Concentration 
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shows  quite  a  different  behavior  (Fig.  2),  with  ethylene  oxide  rather 
intermediate  between  these  examples.  The  decreases  in  pressure  seen 

fupp™*^^  exhaustion  of  the  oxygen 

Figures  3  and  4  illustrate  the  rather  small  effects  of  charge-to- 
mass  ratio  on  the  overpressure.  Figure  3,  for  Comp.  B-3,  is  typical  of 
the  oxygenated  materials,  with  decreasing  pressure  at  low  total  concen¬ 
tration  and  increasing  at  high.  The  opposite  effect  is  seen  in  Figure 
4,  for  THDC.^  At  low  concentrations  (small  values  of  C  +  M)  the  systems 
are  oxygen-rich,  virtually  complete  combustion  occurring  in  all  cases. 
The  heating  value  (AU  of  combustion)  per  unit  mass  is  about  50%  greater 
for  a  hydrocarbon  than  for  magnesium  (due  to  the  lower  atomic  mass  of 
carbon);  hence  increasing  C/M  (i.e.,  replacing  magnesium  fay  fuel)  re¬ 
sults  in  increased  temperature,  and  thus  pressure,  with  the  hydrocarbon. 
But  with  the  explosive  fuels,  the  heating  value  is  less  than  magnesium 
per  unit  mass  (a  result  of  the  nitrogen  and  oxygen  content)  so  that  a 
decrease  is  seen  in  this  same  region.  For  ethylene  oxide,  which  has 
about  the  same  heating  value  as  magnesium,  there  is  very  little  effect 
of  C/M  on  overpressure  at  low  concentrations. 

At  high  concentrations  the  systems  are  all  oxygen— deficient .  Re¬ 
placing  magnesium  by  hydrocarbon  (THDC)  serves  to  exacerbate  the  oxygen 
deficiency,  whereas  the  oxygen  content  of  any  of  the  explosive  fuels 
chosen  will  improve  the  oxygen  balance  and  thus  Increase  the  pressure. 

In  Figure  5  is  shown  the  effect  of  the  different  fuels,  all  at 
C/M  =  1.0  and  for  a  variety  of  total  concentrations.  The  ordinate  of 
this  plot  is  a  "relative  oxygen  content,"  computed  for  each  pure  fuel 
as  the  number  of  oxygen  atoms  divided  by  the  sum  of  the  number  of  carbon 
atoms  plus  half  the  number  of  hydrogen  atoms.  For  the  most  part,  there 
are  no  profound  effects  to  be  seen;  there  is,  however,  a  broad  maximum 
in  pressure  at  HMX,  which  is  fuel  just  oxygen-balanced  to  CO  plus  H2O. 

At  the  higher  concentrations,  ethylene  oxide  shows  a  larger  overpressure 
than  might  be  expected.  This  is  a  result  of  its  high  hydrogen  content, 
leading  at  these  conditions  to  a  large  partial  pressure  of  hydrogen  gas. 
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FIGURE  3.  Comp.  B*3:  Overpressure  vs  Charge-to-Mass  Ratio  at  Various  Total  Concentrations  {in  Kg/m^).' 
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Appendix  A 

EQUILIBRIUM  AND  ATOM-BALANCE  RELATIONS 


The  thermodynamic  equilibrium  constant  for  a  reaction  involving 
gases  is  expressed  in  terms  of  the  partial  pressures  (for  ideal  gases- 
fugacities  in  general)  of  the  reactant  and  product  species.  This  is 
to  say  that  the  standard  state  for  the  gas  is  at  unit  pressure  (fugac- 
ity).  Such  equilibrium  constants,  denoted  Kp,  are  those  usually  found 
tabulated  in  the  literature. 

For  calculations  in  constant-volume  systems,  however,  it  is  more 
convenient  to  use  the  equilibrium  constant  expressed  in  mole  numbers, 
Kj^.  Writing  a  chemical  equation  as  ’ 


i 

where  Aj^  are  the  formulas  of  the  chemical  species  as  the  stoichio¬ 
metric  coefficients,  negative  for  reactants  and  positive  for  products. 


.n 


Vi 


Kp  1  Pi  (assuming  all  substances  are  ideal  gases) , 


.n 


n  ^  and  it  follows  from  the  ideal  gas  law 


that 


V 

RT 


~\  Av 


where  Av  is  the  change  in  the  number  of  moles  of  gas.  The  inclusion  of 
pure  condensed  phases  causes  no  complication;  each,  when  present,  is  in 
Its  standard  state  (unit  activity)  and  it  is  omitted  from  the  equili¬ 
brium  expression  and,  likewise,  not  counted  in  the  evaluation  of  Av. 

The  data  of  Table  3  yields  ^  values  (taken  from  References  1  and 
2;.  Save  for  the  vapor  pressure^  of  condensed  MgO  (the  last  two  given) 
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these  are  all  equilibrium  constants  of  formation;  i.e.,  for  the  syn¬ 
thesis  of  one  mole  of  the  indicated  substance  from  the  elements  in 
their  reference  states:  H2(g),  graphite,  02(g),  N2(g),  Mg(g).  The 
selection  of  magnesium  vapor  as  the  reference  state  for  magnesium,  even 
at  temperatures  below  the  normal  boiling  temperature,  1378  K,  was  made 
for  reasons  of  convenience  in  fitting  the  data  and  leads  to  no  incon¬ 
sistencies  in  results. 

Table  4  lists  values  for  the  logj^Q  Kp  of  the  dissociation  of 
Mg3N2(s)  into  gaseous  elements: 

Mg3N2(s)  =  3Mg(g)  +  N2(s) 

that  is. 


Data  above  1300  K  are  taken  directly  from  Reference  3;  those  below  that 
temperature  include  also  the  vapor  pressure  of  solid  or  liquid  magnesium 
from  the  same  source. 

Formulas  enclosed  in  brackets  represent  mole  numbers;  the  symbol 
a^^  represents  the  number  of  moles  of  atoms  of  element  X  in  the  system. 
The  six  atom  balance  equations  are: 

[Ar]  =  0,4036 

[CO]  +  [CO2]  =  [C(s)]  = 

[H]  +  [OH]  +  2[H20]  +  2[H2]  =  an 

-[NO]  +  [N2]  =  ajj 

[OH]  +  [H2O]  -b  X^MgO  -I-  [NO]  +  [0]  -b  2[02]  =  a^ 

[Mg]  -b  ^MgO  =  ajig, 

where  ^MgO  =  [Mg0(s  or  1)  ]  -b  [MgO(g)]. 

In  actual  practice  the  oxygen  and  magnesium  balances  are  often 
combined : 


ao  -  aj^g  [OH]  -b  [H2O]  +  [NO]  -b  [0]  -b  2  [03]  -  [Mg] 

thus  eliminating  the  need  to  compute  the  mole  numbers  of  condensed  MgO. 
Moreover,  in  the  computational  routines,  what  we  may  refer  to  as 
"adjusted"  stoichiometric  conditions  are  used.  For  example, 

[H2O]  -b  [H2]  =  1/2  ("ajj  -  [H]  -  [OH] 
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making  use  of  previously  computed  or  estimated  values  of  [H]  and  [OH] 
Similar  adjusted  oxygen— magnesium  balance  relations  are  used  in  a 
variety  of  forms. 
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Appendix  B 

EQUILIBRIUM  COMPUTATIONS 


The  following  is  an  outline  of  the  computational  algorithm  re¬ 
ferred  to  in  the  body  of  this  report  in  the  section  "Equilibrium  Con¬ 
siderations".  The  notation  is  that  described  in  Appendix  A. 

[Ar]  is  given  by  its  value  in  air,  0.4036  moles/m^  and  [MgO(g)] 
is  computed  from  the  vapor  pressure  of  the  appropriate  condensed  phase. 
The  mole  numbers  of  the  condensed  phases  themselves  are  computed  after 
all  other  computations,  making  use  of  the  atom  balance  equation: 

[Mg0(s  or  1)]  =  aj^  -  [Mg]  -  [MgO(g)] 

When  present, 

[C(s)]  =  ag  -  [CO]  -  [CO2I 

In  each  iteration,  after  the  principal  computations  have  been 
done  the  "high-temperature  species"  H,  OH,  0,  and  NO  are  found  from 
the  appropriate  equilibrium  constants  of  formation  (see  Table  1).  [N2] 

is  then  adjusted  for  the  loss  of  N  atoms  into  NO. 

There  remain  6  gaseous  species  (referred  to  in  the  following  as 
the  "major"  species)  to  consider:  H2,  H20>  ^0,  CO2,  O2,  Mg. 

Case  I:  Oxygen  is  in  excess  over  both  C  (to  CO)  and  Mg;  that  is, 
^0~^C“^Mg  positive.  [02^0  assumed  equal  to  this  excess;  the 
present  value  of  [O2]  is  used  to  solve  the  pair  of  simultaneous  equa¬ 
tions  involving  carbon-atom  balance  (CO  and  CO2  only,  since  solid  car¬ 
bon  will  not  be  present  in  this  case)  and  the  CO  -  CO2  equilibrium  ex¬ 
pression: 

[C02]/[CO][02]^''^  =  KcojKco 


In  like  manner  [Hj]  and  [H2O]  are  computed  from  the  adjusted  hydrogen 
balance  (see  Appendix  A)  and  the  formation  equilibrium  constant  of  H2O. 
[Mg]  is  computed  from  the  known  [MgO]  and  the  present  approximation  to 
[O2].  Following  computation  of  the  high-temperature  species,  a  new 
[O2]  is  found  from  the  atom  balance  in  oxygen;  this  is  averaged  with 
the  preceding  value  to  compensate  for  overcorrection.  Convergence  is 
tested  on  an  adjusted  oxygen-balance  condition. 
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Case  II :  Oxygen  is  insufficient  for  both  carbon  and  magnesium;  no 
solid  carbon  is  present. 

The  initial  [O2]  is  zero;  successive  approximations  are  generated 
at  the  end  of  each  iteration.  The  remaining  major  species  are  solved 
through  use  of  three  atom-balance  conditions:  C,  H  (adjusted)  and  0 
(adjusted);  and  two  equilibrium  relations: 

[C02][H2]/([C0HH20]  =  Kco2 /(Kc0%20) 

and 

[C02][Mg]/([C0]  [MgO]  =  Kco^ /(Kco^Mgc) 

recalling  that,  in  the  last,  [MgO]  is  a  known  quantity.  These  five 
simultaneous  equations  lead  to  a  cubic  equation  in  [002],  which  is 
solved  by  the  Newton-Raphson  approximation  method.  The  correct  root 
will  lead  to  positive  values,  on  back  substitution  into  the  system  of 
simultaneous  equations,  for  all  mole  numbers.  [O2]  is  computed  from 
[Mg]  and  the  formation  equilibrium  constant  of  MgO.  After  checking  an 
adjusted  oxygen  balance,  this  new  [O2],  averaged  with  the  previous  one, 
is  used  to  initiate  the  next  approximation. 

After  convergence,  the  presence  or  absence  of  solid  carbon  is 
checked  by  evaluating  [CO2]  /  ( [02]K<;;02)  •  ^hls  is  greater  than  unity, 
solid  carbon  is  present  and  the  logic  switches  to  Case  III  below. 

Case  III:  Solid  carbon  is  present. 

if  3-Mg  “  ^0  positive,  this  difference  is  used  as  the  initial 
approximation  for  [Mg].  [O2]  is  computed  from  the  known  [MgO],  Kjjgg 

and  the  present  approximation  for  [Mg].  From  [O2]  can  be  computed  qcO] 
and  [CO2],  using  the  appropriate  formation  equilibrium  constants,  and 
also  [H2O]  and  [H2]  by  Including  the  adjusted  hydrogen-atom  balance. 

A  new  [Mg]  is  found  from  the  combined  magnesium-oxygen  balance  relation 
for  the  next  iteration. 

If  af^  -  Sq  is  negative,  it  is  necessary  to  consider  the  equilibrium 
C0(g)  +  Mg(g)  =  C(s)  +  MgO(g) 


that  is, 

[CO]  [Mg]  =  [MgO]  Kco  %go 
and  an  adjusted  oxygen-magnesium  balance: 

w  =  ag  -ajyjg  -  [OH]  -  [0]  -  [NO]  -  2[02]  -  2[C02]  -  [H2O]  =  [CO]  -  [Mg] 
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The  solution  to  the  quadratic  equation  resulting  from  this  simul¬ 
taneous  pair  yields  [Mg],  from  which  all  the  remaining  mole  numbers 
are  computed  as  in  the  first  statement  in  this  case.  A  corrected  value 
of  w  is  used  for  the  iteration. 


At  the  conclusion  of  computation,  a  check  was  made  for  the  pos¬ 
sible  formation  of  crystalline  magnesium  nitride.  The  product  PMg_PN2 
was  computed;  if  the  computed  product  is  greater  than  the  equilibrium 
value  (Table  4)  the  indication  is  that  Mg3N2  Indeed  forms  and  the  pre¬ 
vious  results  for  that  point  are  discarded 
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Appendix  C 

CONDENSED  STORAGE  OF  INTERNAL-ENERGY  AND 
EQUILIBRIUM-CONSTANT  PARAMETERS 


Memory  space  in  the  WANG  computer  is  severely  limited  for  a  pro¬ 
gram  of  this  complexity:  there  are  248  memory  registers  available  for 
storage  of  data  or  program,  each  two  program  steps  requiring  16  memory 
registers.  To  help  overcome  this  limitation  a  condensed  form  of  stor¬ 
age  for  the  thermochemical -thermodynamic  parameters  was  used.  Thus 
two  registers  were  used  to  store  the  five  internal-energy  parameters 
of  Table  2  (with  sign)  for  each  substance;  and  a  single  register  per 
substance  sufficed  for  the  three  equilibrium  parameters  of  Table  3.  To 
give  one  example  only,  the  constants  A  and  C  of  the  internal-energy 
equation  (Uj  =  A  +  Bt  +  Ct^  +  Dt^  +  E  In  t)  are  stored  in  the  form 
a.c.,  where  a  =  A  (the  sign  of  A  being  given  by  the  sign  of  the  stored 
number)  and  O.c  =  10  =(C  +  10^),  allowing  thus  for  negative  as  well  as 
positive  values  of  C.  The  only  limitation  to  this  method  lies  in  the 
number  of  significant  digits:  for  the  entire  contents  of  the  storage 
register  twelve  digits  are  stored,  not  counting  sign,  decimal  point, 
or  power  of  ten. 
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